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The case against the progenitor's carbon-to-oxygen ratio as a 
source of peak luminosity variations in Type la supernovae 
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Abstract. One of the major challenges for theoretical modeling of Type la supernova explosions is to explain the diversity of 
these events and the empirically established correlation between their peak luminosity and light curve shape. In the framework 
of the so-called Chandrasekhar mass models, the progenitors carbon-to-oxygen ratio has been suggested to be a principal source 
of peak luminosity variations due to a variation in the production of radioactive '^Ni during the explosion. The underlying idea 
is that an enhanced carbon mass fraction should result in a more vigorous explosion since here the energy release from nuclear 
reactions is increased. It was suspected that this would produce a higher amount of '*Ni in the ejecta. In this letter we describe 
a mechanism resulting from an interplay between nucleosynthesis and turbulent flame evolution which counteracts such an 
effect. Based on three-dimensional simulations we argue that it is nearly balanced and only minor differences in the amount of 
synthesized ^*Ni with varying carbon mass fraction in the progenitor can be expected. Therefore this progenitor parameter is 
unlikely to account for the observed variations in Type la supernova luminosity. We discuss possible effects on the calibration 
of cosmological measurements. 



1. Introduction 

The question of the origin of the observed diversity of Type 
la supernovae (SNe la) has been a long standing problem. 
Although these astrophysical events possess a surprising ho- 
mogeneity of features, which led to the application of SNe la 
as distance indicators in cosmology, variations of photometric 
and spectroscopic data have been reported. Nevertheless, cor- 
rections of the luminosity based on an empirical relati on be- 
tween peak luminosity and light curve shape ( Phillips^ <1993l) 
made SNe la a major tool of obse rvational cosmology (for a 
recent review see lLeibundgu t ^OOO). A theoretical understand- 
ing of this correlation is still lacking. Thus the challenge for 
astrophysical theory is to construct a SN la explosion model 
that is robust and — as a first step — able to explain the observed 
variations. 

It is generally agreed on the fact that SNe la originate from 
thermonuclear explosions of white dwarf (WD) stars, although 
many scenarios have been proposed for the particular realiza- 
tion (see Hillebrandt & Niemever (2000) for a recent review). 
The currently favored model is a binary system in which a WD 
composed of carbon and oxygen accretes matter from a non- 
degenerate companion until it reaches the Chandrasekhar mass. 
At this point, thermonuclear burning at the center of the star de- 
velops into a flame which propagates outward. From the view- 
point of hydrodynamics, two modes of flame propagation are 
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possible: a subsonic deflagration, in which the reaction is me- 
diated by thermal conduction of the d egenerate elect rons, and a 
supersonic shock-induced detonation. lArnettI ( Il969l) noted that 
to produce the observed intermediate mass elements a prompt 
detonation can be excluded. Therefore the flame starts out in 
the deflagration mode and may or may not develop into a deto- 
nation. The key feature of this model is that buoyancy-induced 
instabilities lead to a turbulization of the flame front in the de- 
flagration mode. Only due to this effect the flame is accelerated 
sufficiently to power SN la explosions. 

In the first stages the prevailing densities and temperatures 
are high enough that the reaction terminates in nuclear statis- 
tical equilibrium (NSE), consisting mainly of iron group el- 
ements. As the WD expands, density and temperature drop 
and the thermonuclear reactions produce intermediate mass el- 
ements such as Si, S, and Ca. 

The light curve of SNe la is powered by the radioactive de- 
cay of ^^Ni and ^*Co. The peak luminosity is a measure of th e 
^*'Ni produced by the explosion ("Arnett's law". lArnet3ll982h . 
Parameters that have frequently been suggested to affect the 
amount of ^*Ni synthesized are the carbon-to-oxygen (C/O) ra- 
tio of the WD, its central density prior to ignition, and its metal- 
licity. Of course, other parameters like rotation and asphericity 
of the explosion may play a role, too. A thorough investigation 
of the impact of these parameters is mandatory for the valida- 
tion of cosmological measurements, since here possible evolu- 
tion effects with cosmic age are critical. 
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Fig. 1. Temporal evolution of the burning front. 



The objective of the present study is to explore the effect 
of the progenitor's C/O ratio on the supernova explosion by 
means of three-dimensional hydrodynamical models. Since we 
are mainly interested in the explosion energy and the amount 
of synthesized ^^Ni, it is justified to focus on the deflagration 
stage, leaving aside a possible delayed detonation. Only in this 
first stage of the SN la explosion the prevailing densities and 
temperatures are sufficient for burning to iron group elements, 
and therefore the main part of the ^^Ni and probably also of the 
explosion energy is produced here. 

2. The numerical model 

The numerical model applied to simulate the thermonu- 

clear explosion is the s ame as used for several studies by 
^ 1 1 

Reinecke et a l. ( 2002a b). Therefore we will be short in the de- 
scription of the model and only mention the basic aspects here. 

The vast range of involved length scales makes fully re- 
solved SN la explosion simulations impossible. Therefore we 
describe the flame propag ation bv applying the level set method 
jOsher & SethianI Il988l) . Neither the internal flame structure 
nor its wrinkling on small scales are resolved but the flame 
is rather modeled as a discontinuity separating fuel and ashes. 
This discontinuity represents the mean position of the flame. 
To track its propagation it is associated with the zero level set 
of a scalar fi eld G which is evolve d according to the scheme 
described bv lReineckeet al.l l ll999h . In this scheme the effec- 
tive flame velocity must be provided and for this we employ 
the sub grid scale model proposed by Niemever & Hillebrandt 
( 1 1995b which describes the effects of the turbulent motions on 
unresolved scales. The hydrodynamics is modeled based on the 
PROMETHEUS implement ation (Frvxell et al.. 1989) o f the 
piecewise parabohc method jColella & Woodwardlll984) . 

Due to the restricted computational resources only a very 
simplified description of the nucleosynthesis is possible con- 
current with the explosion simulatio n. We follow the approach 
suggested bv'Rein ecke et al.lll200 2a^■ who include five species, 
namely a-particles, '^C, '^O, ^"^Mg as a representative of in- 
termediate mass elements and ^^Ni as a representative of iron 
group nuclei (denoted as "Ni" in the following to avoid confu- 
sion with the particular isotope as part of the ejecta). The fuel 
is assumed to consist of a mixture of carbon and oxygen. At 
the initially high densities burning proceeds to NSE composed 
of a-particles and "Ni". Depending on temperature and den- 



Table 1. Characteristics of the explosion models. 





[10^»erg] 


M("Ni") [Mo] 


max(M^) [Mo] 


0.30 


8.85 


0.5178 


0.04579 


0.46 


9.46 


0.5165 


0.05177 


0.62 


9.97 


0.5104 


0.05636 



sity in the ashes, the amount of a-particles and nickel changes. 
Once the fuel density drops below 5.25 x lO^gcm^^* due to 
the expansion of the WD, burning is assumed to terminate at 
intermediate mass elements and below 1 x 10^ g cm"^ burning 
becomes very slow and is not followed anymore. 

In this letter we will present three simulations that span one 
octant of our coordinate system and assume mirror symmetry 
to the other octants. The simulations were set up on a cartesian 
computational grid that was equally spaced in the inner regions. 
In order to be able to follow the explosion for a longer period of 
time in the outer parts the width of the grid cells was expanded 
exponentially. 

The resolution of the runs was rather low — the computa- 
tional domain was divided in [256] grid cells corresponding 
to a central grid resolution of 10^ cm. In each direction the grid 
length in the outer 35 zones was increased subsequently by a 
factor of 1 . 15. The burning was centrally ignited and the spher- 
ical initial flame structure was perturbed with three toroidal 
rings per octant (see also Fig.Q. 

3. Results of three-dimensional explosion models 

Snapshots from the explosion simulation forX('^C) = 0.46 are 
given in Fig.[n Here the position of the flame front is rendered 
as represented by the zero level set of the scalar field G divid- 
ing fuel and ash. The spatial extent of the burnt region can be 
inferred from the plotted grid which also visualizes our setup 
with uniform grid cefls in the inner region and an exponential 
growth of the grid spacing further out. The flame as initialized 
in our setup is shown in the left snapshot of Fig.^ In the sub- 
sequent images the growth of instabilities and an increasing 
wrinkling of the flame front is visible. 

To test the effect of a varying carbon mass fraction of the 
WD, we performed three simulations with values of X('^C) 
given in Tabled Figure |2l shows snapshots from models with 
different X('^C) after 1.0 s. It is essential for the following 
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XC^C) = 0.62 




Fig. 2. Flame front att - 1.0 s for models with different carbon 
mass fraction of the progenitor material. 



argumentation to note that for all three models the extent of 
the burnt region and also the flame morphology are surpris- 
ingly similar. The latter i s in agre ement with the findings of 
ICiamezo et alJ ll2003i) . but iKhokhlovl ll2000 ) claims that a de- 
creasing X('^C) would result in weaker explosions because of 
a delay of the development of the Rayleigh-Taylor instability, 
which seems to be only a minor effect in our simulations. 

The energy release of the different models are summarized 
in Tabled Obviously, a higher carbon mass fraction leads to a 
higher total energy production for fixed central densities. This 
trend is not surprising. For larger X('^C) the laminar burn- 
ing velocity increases (Timmes & Wooslev, 1992). This effect, 
however, will be important only in the earliest stages of the ex- 
plosion, since flame propagation becomes determined by turbu- 
lence shortly after ignition and can therefore only account for 
minor changes. More significant, a higher carbon mass fraction 
will increase the total energy generation for the simple reason 
that the binding energy of '^C is lower than that of "'O so that 
it releases more energy by fusion to "Ni". 

Therefore the very similar flame evolution of all models 
is not expected in a simple picture. That this is nevertheless 
the case is corroborated by the following arguments. As men- 
tioned in Sect. ^ the flame evolution in the deflagration mode 
is dominated by turbulization which results from Rayleigh- 
Taylor-like instabilities. An important parameter for this is the 
effective gravity experienced by the flame. It accelerates the 
development of flame structures resulting from the non-linear 
stage of the Rayleigh-Taylor instability. At early stages the en- 
ergy fed into the turbulent cascade from the interfaces of large 
scale buoyant burning bubbles will be enhanced by a higher g. 
Thus, the turbulent cascade will develop more quickly and the 
turbulence-induced boost of the effective flame propagation ve- 
locity sets in earlier 

Our models, however, do not show such a behavior In 
Fig. |3l the temporal evolution of the mean g experienced by 
the flame is plotted. We note only small differences for varying 
X('^C). This results in a very similar evolution of the turbulent 
energy in our models. 

The mechanism that is responsible for the reduced effect 
of the varying carbon mass fraction on the flame propagation 
becomes plausible from the evolution of the chemical composi- 
tion of our models (see Fig.l^. Though the temporal evolution 




Fig. 3. Mean gravitational acceleration ^mean experienced by 
the flame front for models with different carbon mass fraction 
of the progenitor material. 




Fig. 4. Temporal evolution of the chemical composition in 
models with different carbon mass fraction of the progenitor 
material. 



of the mass fractions of carbon and oxygen as well as nickel 
show little differences, the variations of the mass fractions of 
a-particles and intermediate mass elements is significant. 

As pointed out by Reinecke et al. (_2Q02a,) . the abundance 
of ff-particles is important for the explosion dynamics. At high 
densities and temperatures the flame converts the C/O fuel to 
NSE which in our models is represented by a mixture of nickel 
(representing the iron group nuclei) and a-particles. The a- 
particles are produced at high temperatures in the ashes and 
have two effects. First, the binding energy of the ashes is lower 
in case of higher X{a), so less energy is released. Second, the 
number of particles per unit mass of the ashes increases and this 
decreases the temperature. Both effects result in a lower tem- 
perature and a higher density of the burning products which 
delays the expansion of the WD and decreases the buoyancy of 
the burnt regions. Thus the flame acceleration is lower and the 
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burning is suppressed. With further expansion of the WD, the 
a-particles are converted to nickel. 

In our simulations a-particles are present between 0.2 s and 
0.9 s. The maximal X{a) is given in Tabled Obviously, higher 
carbon mass fraction in the fuel gives rise to a larger fraction 
of a-particles in the ashes. Therefore the evolution of the flame 
front is more delayed and consequently the explosion dynam- 
ics of the models in the first stage when the iron group elements 
are synthesized is comparable. The result of this effect is a pro- 
duction of similar amounts of iron group elements. The energy 
stored in the a-particles is released at later times and the total 
energy released in the explosion models varies for about 12% 
while the mass of produced nickel differs only about 1 .4% (cf. 
TableHJ. 

Note that the explosion energies and the masses of syn- 
thesized iron group elements in all our models are rather on 
the low side to explain a prototype SNe la. This is due to the 
low numerical resolution of our mode ls. Although the y are ex- 
pected to be numerically converged jReinecke et all Eo02c ). 
the resolution does not allow to apply multi-spot ignition sce- 
narios which have been shown to produce more vigorous ex- 
plosions (.Reinecke et al.^.2002h) . However, the trends inferred 
from our models are expected to be robust. 

4. Conclusions 

In this letter we described a mechanism that is responsible for 
the somewhat surprising fact that although the progenitor's C/O 
ratio affects the energy release of SN la explosions, it has lit- 
tle effect on the peak luminosity determined by the ^^Ni mass. 
This is found with help of three-dimensional simulations and 
dis agrees with the result s reported from one-dimensional mod- 
els. TlH^c^^OCH^ find a 14% decrease in when re- 
ducing the C/O ratio from 1/1 to 2/3; the fact that they consid- 
ered a delayed detonation model should not greatly affect the 
comparability for reasons given in Sect. ^ The "working hy- 
pothesis" of ar i increased ^''Ni pro duction with higher X('^C) 
established bv lUmeda et alJ ( 1 19991) cannot be confirmed. The 
reason for this discrepancy is a complex interplay between nu- 
cleosynthesis and nonlinear flame evolution which can be mod- 
eled appropriately only in three dimensions. The models pre- 
sented here are based on a rather simplistic description of the 
nuclear reactions, which will be improved in forthcoming mod- 
els. Nevertheless, the features that are important for the explo- 
sion dynamics are taken into account and trends can be revealed 
with this approach. 

An interesting question is how our results affect the peak 
luminosity-light curve shape relation of SNe la. Since no light 
curves were calculated from our models we c an only specu- 
late on the effects based on th e trends found bv lArnettllll982h 
and Pinto & Eastman by means of analytic studies. 

According to 'Arnett's law" the peak luminosity should re- 
flect the amount of radioactive ^^Ni synthesized in the explo- 
sion. Since in our models the evolution of the densities of the 
ashes are very similar, the little variation in the iron group el- 
ements will result in little variation of the ^^Ni mass. Thus the 
peak luminosity should be roughly the same for our three mod- 
els. The explosion energy, however, varies significantly and 



a greater expansion at late stages of the explosion associated 
with a larger C/O ratio would result in a more rapid decline 
of the light curve but also in a somewhat higher luminosity. 
This is in contradiction to the empirically established relation 
that for brighter SNe la the light curve declines slower, which 
has been extensively applied to calibrate cosmological distance 
measurem ents. 

Timm eTet alJ ( l2003h analytically found a linear depen- 
dence of the amount of produced ^^Ni on the progenitor's 
metallicity. This in turn is, however, not likely to affect the 
explosion energy and dynamics significantly. Of course, from 
the viewpoint of stellar evolution there is a interrelation be- 
tween the progenitor's metallicity and other parameters, in- 
cluding its carbon mass fraction. Therefore it is plausible that 
the "luminosity-width relation" is caused by a combination 
of different parameters and that a direct relation between the 
^*Ni mass and the peak luminosity might be oversimplified. 
It seems well possible that the distribution of the ^^Ni due to 
three-dimensional effects plays an important role here. In a 
subsequent publication we will report on a systematic parame- 
ter study addressing these questions. Ultimately, SN la explo- 
sion models will have to be coupled to realistic stellar evolution 
of the progenitor system and the calculation of synthetic light 
curves becomes mandatory. 
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